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Abstract 

The description of the universe evolving in time according to general relativ¬ 
ity is given in comparison with the quantum description of the same universe 
in terms of semiclassical wave functions. The spacetime geometry is deter¬ 
mined by the Robertson-Walker metric. It is shown that the main equation of 
the quantum geometrodynamics is reduced to the non-linear Hamilton-Jacobi 
equation. Its non-linearity is caused by a new source of the gravitational field, 
which has a purely quantum dynamical nature, and is additional to ordinary 
matter sources. In the semiclassical approximation, the non-linear equation of 
motion is linearized and reduces to the Friedmann equation with the additional 
quantum source of gravity (or anti-gravity) in the form of the stiff Zel’dovich 
matter. The semiclassical wave functions of the universe, in which different 
types of matter-energies dominate, are obtained. As examples, the cases of the 
domination of radiation, barotropic fluid, or new quantum matter-energy are 
discussed. The probability of the transition from the quantum state, where 
radiation dominates into the state, in which barotropic fluid in the form of dust 
is dominant, is calculated. This probability has the same order of magnitude 
as the matter density contrast in the era of matter-radiation equality. 

PACS numbers: 98.80.Qc, 98.80.Cq, 95.35.-|-d, 95.36.-|-x 

1 Introduction 

The quantum theory predicts that the wave properties should be observed in the 
matter system which this theory describes. The main object of the theory is the 
state vector (wave function) T. The state vector satisfies the differential equations 
defined in some configuration space If. In the general case, the state vector T is 
a complex-valued function, and without loss of generality it can always be written 
as Tq, = Aa 6*“^“, where Aa and ^pa are real functions of generalized variables in 
q, a is a set of quantum numbers which characterize the state of the system with 
the state vector Tq. In the region of kl, where the phase pa varies by a large 
amount on small scales, the system under investigation can be considered as an 
almost classical system, in a sense that its wave properties are inessential and can 
be ignored when calculating parameters of the system [T]. Nevertheless, an almost 
classical system still has the wave properties. They can give probabilistic character 
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to parameters of the system which have no analogs in classical theory. So, the overlap 
integral, ('ho'l'I'a) = f dVlA*^iAa at a' A can be nonzero because of the 

contribution of a subregion in Q, where the difference ipa — 'Pa' of two large phases is 
small. As a result, the spontaneous transition (or transition under the action of an 
instantaneous perturbation) a —)• a' with the change of physical state of the system 
becomes possible. Classical and quantum descriptions of physical properties of the 
same system appear here as complementary without contradiction to each other. 

In the present paper, we provide the comparative description of the evolution of 
the universe as a whole in terms of classical and quantum theories. As is well-known, 
the universe is subject to classical laws of general relativity on large spacetime scales, 
whereas on small, comparable with Planck scales, it should be described from the 
quantum-theoretical perspective. The questions can be posed whether the universe 
preserves the wave properties during its subsequent evolution and whether these 
properties can be discovered. 

As a working model of spacetime geometry we choose an example of a minisuper¬ 
space model with the maximally symmetric geometry described by the Robertson- 
Walker metric. Despite its relative simplicity, this model may give insights into 
understanding of a complete theory of quantum cosmology. In Sect. 2, the quantum 
constraint equations imposed on the state vector of the universe are given. These 
equations are formulated in the representation of the generalized field variables, 
such as the cosmic scale factor and the uniform scalar field. The scalar field is 
described by some Hermitian Hamiltonian. Its mean values with respect to proper 
state vectors determine the proper energy of matter in the form of barotropic fluid 
contained in the comoving volume (Sect. 3). As is shown in Sect. 4, the main equa¬ 
tion of the theory can be rewritten as the nonlinear Hamilton-Jacobi equation. Its 
nonlinear part is caused by a new source of the gravitational held, which has a 
purely quantum dynamical nature, and is additional to ordinary matter sources. In 
Sect. 5, the classical description of the universe evolving in time according to power 
and exponential laws are given in comparison with the quantum description of the 
same universe in the language of wave functions in the semiclassical approxima¬ 
tion. The semiclassical wave functions of the universe, in which different types of 
matter-energies dominate, are obtained. As examples, the cases of the domination of 
radiation or barotropic fluid are discussed. The case of domination of new quantum 
matter-energy is special. It is shown that its energy density is negative, while its 
equation of state coincides with the equation of state of the stiff Zel’dovich matter. 
Such an energy density dominates in the sub-Planck region. Here the wave function 
is constant, and the semiclassical equation of motion has an allowed trajectory in 
imaginary time. The fact that the wave function is non-vanishing near the initial 
singularity point means that in this region there is some source which provides the 
origin of the universe with a finite nucleation rate [21 [3] (cf. Refs. [H O [6]). In 
Sect. 6, as a supplement, the transition probability of the universe from the state, 
where radiation dominates, into the state, in which barotropic fluid in the form of 
dust is dominant, is calculated. 

Throughout the paper, unless otherwise specified, the modified Planck system of 
units is used. As a result, all quantities in the equations become dimensionless. The 
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length Ip = ■sj2Ghl (Sttc^) is taken as a unit of length and the pp = 3c^/(SvrGZp) is 
used as a unit of energy density and pressure. The mass-energy is measured in units 
of Planck mass, mpc^ = hc/lp. The proper time r is taken in units of Ip. The time 
parameter (conformal time) T is expressed in radians. The scalar field is taken in 
(pP = -^3c^/(S ttG). Here G is Newton’s gravitational constant. 

2 Quantum constraint equations 

We consider a canonical approach to the quantization of general relativity, in par¬ 
ticular quantum geometrodynamics. It is well-known that approaches of this type 
share a common problem known as the ‘problem of time’ (see, e.g. Refs. Oil]). The 
apparent manifestation of this problem is the absence of explicit time parameter 
in the Wheeler-DeWitt equation considered as the main dynamical equation of the 
theory. 

The way to solve this problem could be rewriting the classical constraint equa¬ 
tions to obtain a Schrodinger-type equation, as a preliminary to quantization. Since 
it was proved that general relativity could not be viewed as a parametrized held 
theory [9], a concept of matter clocks and reference huids was proposed [8]. One 
approach in this direction goes back to DeWitt, who studied a coupling of clocks 
to an elastic media m- It can be shown that perfect huids are a special case of 
DeWitt’s relativistic elastic media and that in the velocity-potential formalism for 
perfect huids the thermasy can be connected to a clock variable with DeWitt-type 
coupling m- Kijowski et al. have proposed a physical interpretation for the ther¬ 
masy as a “material time” relating it to “proper time retardation” due to the chaotic 
motion of the particles of the huid moving chaotically around the how lines [111113]. 
The perfect huid plays a crucial role in the formalism as the “reference huid”, whose 
particles identify spacetime points, and clocks carried with them identify instants 
of time (cf. Ref. [8]). 

We are working in the framework of the velocity-potential version of perfect- 
huid hydrodynamics formulated by Seliger and Whitham m, and generalized by 
Schutz m (see also US]). In this formalism, the perfect huid four-velocity 
is written as a combination of (at least) hve scalar helds (potentials) and their 
gradients. An independent physical interpretation can be given to each potential. 
So, one potential appears to be the thermasy 0 introduced earlier by van Dantzig 
HZ!; another is the entropy s, and the other one is the potential A for the specihc 
free energy and so on. In the general case, all hve potentials are needed in order to 
describe a relativistic perfect huid. In the cosmological model with the Robertson- 
Walker metric and irrotational hows of a perfect huid, the number of potentials 
reduces, and the perfect huid can be described by three potentials: 0, s, and A. 

Let us consider the homogeneous and isotropic cosmological model and write the 
Robertson-Walker line element in the form 

ds'^ = a^[dT^ - dnl], ( 1 ) 

where a is the cosmic scale factor which is a function of time, T is the time variable 
connected with the proper time r by the diherential equation dr = adT, T is the 
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“arc-parameter measure of time”; during the interval dr, a photon moving on a 
hypersphere of radius a(r) covers an arc dT measured in radians [18]. dVL^ is a line 
element on a unit three-sphere. Following the ADM formalism naiso] , one can 
extract the so-called lapse function N, that specifies the time reference scale, from 
the total differential dT: dT = Ndr], where r] is the “arc time” which coincides with 
T for N = 1 (cf. Refs. |18l I21j i. In the general case, the function N plays the role 
of the Lagrange multiplier in the Hamiltonian formalism and it should be taken into 
account in an appropriate way. 

Having in mind mentioned above, we can write down the Hamiltonian for the 
cosmological system (universe) in the form |22l l23l 


H = 


+ 


N 

y 

Ai 


{-T^l - + a^[p^ + p^]] 



Pos 


-|- A2 


® ^0 


( 2 ) 


where vTa, vr©, are the momenta canonically conjugate with the variables a, 0, A, 
is the energy density of matter (the field (/>), p^{po,s) is the energy density of a 
perfect fluid, which defines a material reference frame uni EH, and it is a function 
of the density of the rest mass po and the specific entropy s; the 0 is the thermasy 
which defines the temperature, T = the is the four-velocity; the A 

is the potential for the specific Gibbs free energy T taken with an inverse sign, 
= —X^uU^. The N, Ai, and A 2 are the Lagrange multipliers. 

The Hamiltonian (|2|) is a linear combination of constraints (expressions in braces) 
and thus weakly vanishes, Lf ~ 0. The variations of the Hamiltonian with respect 
to N, Ai, and A 2 give three constraint equations, 

- + a^[p^ + p^] « 0, tt© - ^ a^pos ^0, + \ ~ 0- (3) 

From the conservation of these constraints in time, it follows that the number of 
particles of a perfect fluid in the proper volum^i] and the specific entropy con¬ 
serve: Eq = ^a^po = const, s = const. Taking into account these conservation laws 
and vanishing of the momenta conjugate with the variables po and s, vTp = 0 and 
TTg = 0, one can discard degrees of freedom corresponding to these variables, and 
convert the second-class constraints into first-class constraints in accordance with 
Dirac’s proposal [25] . 

It is convenient to choose the perfect fluid with the density p-y in the form of 
relativistic matter (radiation). Then, in Eq. ([3]) one can put a^^p^ = E = const. The 
matter field with the energy density and pressure can be taken for definiteness 
in the form of a uniform scalar field 4>, 

p<l> = + P4, = ^TTl-V{(j)), ( 4 ) 

where V (</>) is the potential of this field, vr^ is the momentum conjugate with cj). 
After averaging with respect to appropriate quantum states, the scalar field turns 
into the effective matter fluid (see Ref. [23], and below). 

^This volume is equal to 27 r^a®, where a is taken in nnits of length. 
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In quantum theory, first-class constraint equations ([3]) become constraints on the 
state vector 'h [25] and, in this way, define the space of physical states, which can be 
turned into a Hilbert space (cf. Ref. [2B|). Passing from classical variables in Eqs. 
(|2|)-Q to corresponding operators, using the conservation laws, and introducing the 
non-coordinate co-frame 


hdr = sdB — dX, hdy = sdQ + dX, (5) 

where h = is the specihc enthalpy, is the pressure of radiation, and y is a 

supplementary variable, we obtain 


= 0 , dy^ = 0 , ( 6 ) 

(^-dl + a^- 2aH^ - h) ^ = 0, (7) 

where 

( 8 ) 

is the operator of Hamiltonian of the scalar field (j) which is Hermitian, and the 
operator is described by Eq. (|3|) with Equation (|7|) is the Wheeler- 

DeWitt equation of a minisuperspace model for the universe filled with a scalar field 
and radiation, when the state vector 'h does not depend on time. In the approach 
under consideration, the dependence of the state vector 'I' on time is determined 
from Eq. Q. 

The quantum constraints Q and (|7|) can be rewritten in the form of the time- 
dependent Schrodinger-type equation 

— idT^ = (9) 

o 

where 

'H = -dl + a^ - 2aH^. ( 10 ) 

The minus sign before the partial derivative dx is stipulated by the specific character 
of the cosmological problem, namely that the classical momentum conjugate to the 
variable a is defined with the minus sign [271128j (see below). 

The partial solution of Eqs. dH) and ([7|) has a form 

^(T) = e*i^('^"^°)T(ro), (11) 

where the vector T(To) = (a, satisfies the stationary equation 

=E|V’). (12) 


Erom the condition 


0 = ^ / D[a, 0] |T|2 = f D[a, <(>] 


rO -n 




(13) 


where D[a, (j)\ is the measure of integration with respect to the fields a and cj) chosen 
in an appropriate way, it follows that the operator (fTOll is Hermitian: T-L =1-0. 
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3 Barotropic fluid 


The Hamiltonian of matter Hfj, can be diagonalized by means of some state vectors 
{x\uk) in the representation of generalized field variable x = x{a, 4>) with the measure 
of integration D[a, c^] = dadx in Eq. (fT^ . 

Assuming that the states \uk) are orthonormalized, {uk\uk') = dkk', we obtain 
the equation 

{uk\H,p\uk') = Mk{a)5kk', (14) 

which determines the proper energy Mk{a) = ^a^pm of a substance (barotropic 
fluid) in discrete and/or continuous fcth state in the comoving volume with the 
energy density pm = {uk\p<j>\uk) and pressure pm = = Wm{a)pm, where 


Wm{a) 


1 dlnMk{a) 
3 din a 


(15) 


is the equation of state parameter. 

In order to obtain the general form of the proper energy Mfc(a), let us suppose 
that the interaction of the uniform scalar field is described by the potential V{4>) = 
Aqi/)", where is the coupling constant and a takes arbitrary non-negative values, 
a > 0. Then we find 

2 

\ 2+a 3(2-a) 

-f) ( 16 ) 

where is an eigenvalue of the equation (—d^ + — ek)luk} = 0, and x = 

1 

4> is the rescaled matter scalar held. The equation of state parameter 
in such a model does not depend on a and has a simple form 



Wm{a) 


a — 2 
a + 2 


It describes the barotropic fluid in all possible states. In the case of the model cjP, 
the field (p averaged over its quantum states reproduces vacuum (dark energy) in 
the fcth state with the density pm = AoCfc and the function {x\uk) in the plane wave 
form with the wave vector k = — 1. The model describes the strings 

in the fcth state with the energy density Pm = where ^ 0 and \uk) is 

the Airy function. In the model cp"^, the scalar field, after averaging over quantum 
states, turns into dust with the total mass Mk = '/2^{k + ^), where k is the number 
of dust particles, and the density pm = The model cp'^ leads to the relativistic 

matter with the energy density Pm = ^5 where e*, < oo and \uk) has the 

asymptotics in the form of cylindrical function. In the case a = oo, the field (p, 
averaged over the states \uk)-, reduces to the stiff Zel’dovich matter with the density 

_ 2efc 

Pm - 
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4 Non-linear Hamilton-Jacobi equation 


Assuming that the set of vectors \uk) is complete, Ylk \'^k){uk\ = 1, the solution of 
Eq. (|12p can be represented in the form of the superposition of states of the universe 
with the substance in the kth. state in any form described above. We have 

I'fp) = '^\uk){uk\tp), (17) 

k 

where the wave function /(a) = (ukltp) satisfies the equation 

[-dl + a^- 2aM{a)] f = Ef. (18) 

The index k is omitted here and below, since in what follows we consider the universe 
with the proper energy of the substance in a specific kth state, Mk{a) = M{a). 
Because the operator (11011 is Hermitian, it follows that the operator on the left-hand 
side of Eq. (|18p is Hermitian as well. This equation determines the wave function 
corresponding to the particular eigenvalue E. Depending on the form of M{a), the 
constant E can take the values lying in a discrete or continuous spectrum of the 
states of the proper energy of radiation My = with the energy density Pj = ^ 

and pressure = |/3y. Thus, the value of the constant E is determined through 
the quantum numbers enumerating the states of the substance and radiation. 

We look for the solution of Eq. (jl8|l in the form of the wave propagating along 
the a direction 

/(a)=yle*^(“), (19) 

where A is the normalizing constant, and the phase S{a) is a complex function 


S{a) = SR{a) -t- iSi{a) 


( 20 ) 


(with Sr and Sj real). In the general case, the solution of Eq. (fTHI) is the super¬ 
position of the wave function /(a) and its complex conjugate f*{a). Substituting 
Eq. (fT9]l into Eq. (fT8]l . we find that SR{a) satisfies the non-linear equation 

idaSRf + a‘^-2aM{a)-E = Q{a), (21) 


where the function 

3 

4 2daSR 

describes the new source of the gravitational field with the energy density 


PQ 



( 22 ) 


(23) 


which is additional to the ordinary matter (substance and radiation). This source 
has the quantum dynamical nature. It emerges as a result of expansion (or con¬ 
traction) of the universe as a whole. The equation of state of quantum source of 
matter-energy has a form 

PQ = WQ{a)pQ, (24) 
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where pg is the pressure, and 


wgia) 


1 / dlnQ(a)\ 
3 \ din a ) 


(25) 


is the equation of state parameter. The first term in Eq. ()25p takes into account 
the correction for relativity. The additional source of matter-energy has a proper 
energy Mq = \a?PQ = ^ contained in the volume 

The derivative of the imaginary part of the phase S{a) is 


daSi 


I dlSn 

2 daSn 


(26) 


In order to clarify the physical meaning of the quantum correction in Eq. (I2ip . we 
rewrite the source function Q{a) as 


Q{a) = {daSif - dlSi, 


(27) 


where Eq. (I26p was used. It follows from here that if the imaginary part Sj of the 
phase (1211 is a slowly varying function of a, then one can set Q{a) ^ 0. In this case, 
Eq. (I2ip becomes the Hamilton-Jacobi equation for the classical action S^i and the 
wave function (jI9p takes the form 

/(a) = const (28) 


It describes the quantum universe in the semiclassical approximation, when Sji = 
Sci- In classical mechanics, the momentum is equal to the first derivative of the 
action with respect to the generalized coordinate. In general relativity, Hamilton’s 
equations of motion lead to [28] 

daSci = ~^ = ( 29 ) 

where we denote d = ^. Equations ([28P and (l2^ describe the same universe, but 
from the different point of view, namely as the quantum system in the approximation 
Sr = Sci or as the classical object obeying the laws of general relativity. 

To determine the physical meaning of the derivative daSR, we calculate the 
probability flux density for the universe to be a hypersurface with radius a in a 
four-dimensional space. From Eq. m, it follows that the probability flux density 
is described by the expression 

Ja = ^{rdaf-fdan- (30) 

2^ 

Taking into account Eq. (fT^ . we find 

Ja = \f\^daSR, (31) 

where |/p = is the probability density of the universe to have the scale 

factor a. Equation (|3ip shows that the wave function (|19p describes the expansion 




of the universe as a whole with the generalized momentum daS^. The generalized 
action Sr is the solution of the non-linear Hamilton-Jacobi equation (1211) . 

One can make sure that the energy density (1231) is the quantum correction to the 
energy density of the substance and radiation by rewriting Eq. (|2ip in dimensional 
physical units 


where 


{daSnY + 


Pm — 


2 , I Svrc^ ^ ^2 


2G 


M{a) 


1 - 


8ttG 


{Pm + P"! + Pq) 


= 0 , 


_ ^ 


27r2a3’ '^3' a4’ 

are the energy densities of the substance, radiation, and quantum addition measured 
in GeV/cm^. Here a is taken in cm, M{a) in GeV, E in GeV cm (he), Sr in GeV s 
(h), whereas Q is in cm“2 and it has the same form as in Eq. ()22n . Erom Eq. (1321) . 
it follows that the quantum correction is proportional to In the formal limit 
h ^ 0, Eq. (I32p turns into the Hamilton-Jacobi equation of general relativity. A 
more rigorous approach requires the change to dimensionless variables which do 
not contain dimensional fundamental constants G, c, and h. The impact of the 
quantum correction Q{a) on the dynamics of the universe as a whole is determined 
by how quickly the amplitude of the wave function (I19p changes with a. In 

Eq. (I2ip . this impact depends on how its terms behave with increasing (decreasing) 
of o. In the models, in which Eq. (I18p can be integrated exactly, Eq. (1211) also admits 
a solution in an analytical form 


GQ{a) 


(32) 


(33) 


5 Classical-quantum correspondence 

In order to determine the correspondence between classical and quantum description 
of the same physical system, we consider the specific model of the universe. Let us 
assume that the universe evolves in time r with a power-law scale factor 


a = /3T'^, 

(34) 

where a and (3 are some arbitrary constant^. Then the generalized momentum in 

the semiclassical approximation is equal to 


1 2a-l 

OaSR = -a^o-a - , 

(35) 

and the quantum source function (I22p takes the form 


Q[a) = -2 , 

(36) 

where the numerator 


(2a - l)(4a - 1) 

“ 4a2 

(37) 


^Here the constant a differs from the parameter a in Eq. m- We use the same letter to 
emphasize that, in both cases, different types of matter-energy correspond to different values of a. 
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does not depend on a. Universal dependence of the source function (f36|) on the scale 
factor allows one to find the equation of state for quantum matter-energy for any 
values of the parameters a and (5. So, it follows from Eq. (12^ and pUjl that the 
equation of state parameter ([ 2 ^ is wq = 1 , and the energy density decreases with 
the increase of a according to the law 

PO = 5 . (38) 

As was mentioned in Sect. 3, the stiff Zel’dovich matter has such a density. The 
energy density of this quantum matter can be negative (5 < a < ^), positive (a > ^, 
and a < |), or vanish (a = and a = |). 

From the point of view of quantum theory, the expansion of the universe with 
the scale factor (l34|l is described by the semiclassical wave function 


/a (a) 


2a-l 

Aa a 2 “ exp 


, a‘^/3 a 3a-l 1 
~^3a - 1 “ “ J 


at ex 7^ 2, 


(39) 


and ^ 

/i (a) = Ai a2“®4". (40) 

3 3 

The complex conjugate function /^(a) corresponds to the wave propagating towards 
the initial cosmological singularity point, a = 0 , and describes the contracting uni¬ 
verse. 

In the limit of infinitely large values of a, we have 


lim 7 „ = 2, (41) 

Q—>-00 

and the quantum correction in Eq. (|2ip equals 

Q{a) = 4- (42) 

The same additional term is produced by the universe expanding exponentially. 
Really, setting 

a = a(0)e^'^, (43) 

where p-u is some constant (e.g. = ^, where A is the cosmological constant), we 

calculate the momentum of the universe 


daSR = -^a^. (44) 

Substituting (|44p into Eq. (j22p . we obtain the quantum correction in the form (1421) . 
The exponentially expanding universe is described by the semiclassical wave function 

/ 00 (a) = p4oo-e-'^“'. (45) 

a 

At infinity this function vanishes, but it diverges at the point o = 0. However, this 
point is inaccessible, since r > 0 in Eq. (|43p . The function (j45p can be normalized 
to a constant. The probability flux density (j3ip for the wave function (1451) is 

Jeo — V Pv I Aqq I . (46) 
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The approximation (l36|) linearizes Eq. (| 2 T|) . and it can be considered as the 
Hamilton-Jacobi equation for the generalized action Sji with the additional source 
of the gravitational field with the energy density Using Eq. pSjl . this equation 
can be easily reduced to the Friedmann equation for the Hubble expansion rate 
H = ^. In the semiclassical approximation, daSpt = —aa, we have 




2 M(o) E ja 1 

Q I A I a o 


(47) 


In addition to the energy density of the substance (~ M a~^) and radiation (~ a~^), 
this equation contains the energy density of the quantum source (~ a“®). 

Let us find out what restriction on the solutions of Eqs. ()21l) and (1471) is imposed 
by Eq. (1341) . Since it is assumed that the universe evolves according to the law 
(j34p with a given a, then, according to the standard model, it means that the 
approximation of a single component domination in the total energy density of 
matter energy p is used. It has a form [29l [30] 

fora~r“. (48) 

Substitution of Eqs. (1351) and (l36l) into Eq. (1211) gives the condition 


a 


13a — 2aM{a) — E = ^aa 


Let a = then 71 = 0, and Eq. (14^ is equivalent to 


2M(a) 


Then from Eq. (1471) . we obtain 


or H = E 
' 2 a 2 / 2 r’ 


(49) 


(50) 


(51) 


if Eq. (|34l) is used. This equation describes the spatially flat universe, in which 
radiation dominates. Hence, the model (IMl) does not contradict Eq. ([T7|) . Such a 
universe expands with constant momentum 

daSR = (52) 

and, according to Eq. (1221) . an additional quantum source of energy is not generated, 
Q{a) = 0. The quantum properties of such a universe are described by the wave 
function 

fi (a) = Hi = Hi . (53) 

2 2 2 

This function is a solution to Eq. (|18p in which the terms and 2aM (a) are omit¬ 
ted, i.e. the quantum universe is spatially flat and contains nothing but radiation. 
Comparing Eq. ([531) with Eq. ([3^ at a = ^, we find the constant 0 of Eq. (IM|) 

0 = ( 2 VE) . (54) 
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This expression coincides with the one that can be obtained directly from the solu¬ 
tion of Eq. (1471) in the approximation specified above. 

The wave function of the continuous state fi (a) = fi (a; E) (1551) can be normal- 

2 2 

ized to a delta-function 


{fi{E)\fi{E'))=6{E-E'). 
2 2 


(55) 


This condition determines the normalizing constant Ai 

2 

factor), 



(up to inessential phase 


(56) 


According to Eqs. (f5^ and (l54l) . the probability flux density (f 8 T]l for the universe 
with the wave function (I53p and the amplitude from Eq. (I56p does not depend on a 
and E, and it equals 

(57) 

In this case, the conservation law is fulfilled, daJi = 0. The minus sign in Eq. (1571) 

2 

shows that the matter flux is directed away from the observer, i.e. matter objects 
(galaxies) move away from the observer demonstrating the effect of expansion of the 
universe. 

Using the wave packet (proper differential) 


h{a;E) = 
2 


rE+5 

Je-s 


dE'fi{a]E') 


(58) 


with the width 2(5 <C 1, (5 > 0, the wave function fiia]E) can be normalized to 

2 

unity as follows: 


(/i(E)|/i(E)) = l. 
2 2 


(59) 


From Eqs. ([551) and (fSHp . it follows that the integral (l5^ does not depend on 5. The 
smaller is (5, the wave packet ([551) reproduces the wave function (IS5p with greater 
accuracy. The wave function (f53P with the amplitude ([56l) can be interpreted as a 
part of the de Broglie wave propagating along a with the momentum daS^ = —y/E. 

In the case o: = |, we have 72 = ^, and the quantum source function (1221) is 


Q{a) 


5 1 
16 ^' 


(60) 


Quantum source is characterized by the positive energy density pg ([23]). The con¬ 
dition ([3Up is written as 


2M(a) E 5 1 1 _ 4 //? 

16 a® 9 \ a 


and Eq. ((47l) gives the Hubble expansion rate 


= 


4 fp 


9 V a 


or H = 


3t 


( 61 ) 


(62) 
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The latter equation describes the spatially flat universe, where non-relativistic mat¬ 
ter with the mass M{a) = M = const dominates. In this case, the constant /3 
is determined by Eqs. (1621) and (I47|) . in which domination of matter is taken into 
account 



(63) 


The universe expands as a whole with the momentum 


daSR = = -Ha^. 


(64) 


The wave function has a form 

/|(a) = ^exp|-z-VWa3/2| = (65) 

This function is the asymptotics of the Airy function at 2oM S> 1. The Airy function 
is the solution of Eq. (|18l) with M{a) = M = const for the spatially flat universe 
(the term in Eq. (fT8|) should be omitted). In the domain a 3> the Airy 
function describes the continuous state with respect to E and can be normalized to 
a delta-function 5{E — E') [T]. The asymptotics (1651) describes the state with E = 0. 
Therefore, it is convenient to normalize it by the condition 


(/|(M)|/|(M')) =5(M-M'), (66) 

where we denote / 2 (a) = In order to calculate the parameters of the 

3 3 

universe in the state (j65|l . the wave function can be normalized to unity 


(/2(M)|/2(M))=1 

3 3 

instead of (1661) . where 

pivl -\-S 

/ 2 (a;M)= / dM'f 2 {a-,M') 

3 Jm-5 3 


(67) 


( 68 ) 


is the wave packet with the width 25 <C 1, and d > 0. The amplitude A 2 in 

3 

Eq. (IS5|) can be found under the assumption that the probability flux density is 
conserved in the expanding universe and it equals to the probability flux density in 
the radiation-dominated era ([57]) . As a result, we obtain 


1 ^ 


2 

3 


2 


1 

2V^' 


(69) 


The wave function (j65p with the normalizing constant (1691) no longer has the form 
of a part of the de Broglie wave. Its amplitude depends on a, and the phase factor 
contains the momentum (j64|) multiplied by 
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The case a = | in quantum description appears to be special. The wave function 
has the form (00]). The parameter (1371) takes the value 71 = —| and the quantum 
source function ( 1221 ) is 


Q(«)- 4 ^ 2 ’ 


so that the energy density ([23]) is negative 




Then, the condition (I49p is written as 


(71) 


(72) 


2M(a) 


E 

H —7 — 


4a6 


9 I o 


(73) 


Equation (071) together with Eq. (1731) defines the Hubble expansion rate 


= 


1 fP 


9 V a 


or H = 


3t 


(74) 


This equation describes the spatially flat universe, in which the stiff Zel’dovich 
matter with the energy density (1721) dominates. From Eq. (1471) . it follows that for 
such a universe 


= - 


1 

4^’ 


(75) 


because the energy densities of a substance and radiation and the curvature term 
should be neglected. This corresponds to the domain of values a < 1 (sub-Planck 
scales). Equations ((741) and ([75l) impose restriction on allowed values of the param¬ 
eter P 

P^ 1 

y = -4' 

If this condition is satisfied, then it means that in sub-Planck region, where the 
energy density (1721) dominates, the wave function ()40p is either constant 


1 P^ 

fi(a) = Ai at - = i—. 
■’3^' 3 2 3 ' 


or increases linearly with the increase of o, 


1 P^ 

fi{a) = Aia at - = 


(77) 


(78) 


In the case of Eq. ([77p . we have /i(0) = Ai = const, i.e. there is a source at 

3 3 

the point a = 0. It may indicate that the universe can originate from the initial 
cosmological singularity point with the finite nucleation rate P rsj |/i(0)|2. In a 
more rigorous consideration, it appears that such an origin occurs not from the 
point a = 0, but from the whole domain of values of a < where the wave 

function has a form (1771) [21 [3]. In this domain, there exists the classical trajectory 
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in imaginary time t = —it + const, which is the solution of Eq. (|47ll . where the 
energy density of a substance and the curvature term are neglected. Near a ~ 0, 
one has 

/ q \ 1/3 

o=^--iTj (79) 

From this solution, it follows unambiguously the condition on /3 which chooses the 
solution (fTTl) . 

In accordance with Eqs. (fMjl . ([l8l) . and (f63]) . in both cases a = ^ and a = |, 
classical theory predicts the existence of the initial cosmological singularity at the 
point a = 0 in the domain of real values of the scale factor. The solution of Eq. (j75jl 
which takes into account the domination of quantum correction (j72p on scales o < 1 
demonstrates that, according to Eq. ([79]), the initial cosmological singularity lies in 
the non-physical region of imaginary values of a and it is inaccessible from the point 
of view of general relativity. Moreover, if one takes into account quantum effects 
in the region a < 1, it will allow one to revise the properties of the universe on 
sub-Planck scales. Quantum-mechanical description of the universe in semiclassical 
approximation admits the possibility of an origin of the universe from the sub-Planck 
domain. As is shown in Refs. 121 El, the origin of the universe is accompanied by a 
change in space-time topology, so that the geometry conformal to a unit four-sphere 
in a five-dimensional Euclidean flat space changes into the geometry conformal to a 
unit four-hyperboloid embedded in a five-dimensional Lorentz-signatured flat space. 
On the boundary, where these two subregions adjoin each other, there is a jump 
with change of metric signature [HE]. 


6 Example 


Let us calculate the probability of the transition of the universe from the state, 
where radiation dominates, into the state, in which barotropic fluid in the form of 
dust is dominant. This probability is determined by the expression 


K/|l/pP 

w (rad. —dust) =-^^— 




(80) 


Using the wave functions (l53]l and ([65]l . we have 


i(/2|/i)i =i^2ri^ir/(i^,M), 

3 2 3 2 


(81) 


where 


I{E,M) = 


1 


1 


[ da a exp < i^/E ( - _ 

h \ \^a,{E,M) ^ 


(82) 


and 


ac{E,M) 


9 E 
8 M' 


(83) 
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The integration in Eq. (j82ll can be performed analytically. As a result, I{E,M) 
will have the form of the sum of terms containing generalized hyper geometric func¬ 
tions. However, for our illustrative purposes it is sufficient to calculate the integral 
in Eq. (|82l) using the method of approximate calculation of overlap integrals of semi- 
classical wave functions mentioned in the Introduction. In the region a ^ a^, the 
exponential in the integrand oscillates rapidly and its contribution into the integral 
is exponentially small [T]. The main contribution into the integral comes from the 
region near a = Uc, where the exponential is almost unity. In this approximation, 
we have 

Q / p \ 

= 

Then, taking into account Eqs. (1561) . ([59|) . (l67|) . ([691) . and ([8^ . we obtain the fol¬ 
lowing simple expression for the probability (ISOh 


w (rad. —)• dust) = 


3 E 
8 M2’ 


(85) 


We estimate Oc and w using the values Eq = 1.86 x 10^^® and Mq = 0.92 x 10®^. They 
correspond to the modern values of the energy densities of radiation and matter, 
= 2.61 X 10“^° GeV cm“^ and = pcrit = 0.48 x 10“^ GeV cm“^, and the 
Hubble length uq = ^ = 1.37 x 10^® cm taken as a rough estimate of the size of 
the observable universe H. We have 


ac{EQ, Mq) = 2.27 x 10^^(= 1.69 x 10^^ cm), (86) 

wq = 0.83 X 10"^. (87) 

The redshift Zc = 2Alzeq, where Zeq = 3360 is the redshift of matter-radiation 
equality, corresponds to the scale factor Oc ([86]). The probability (l87l) has the same 
order of magnitude as the matter density contrast ~ 10“^ in the era of matter- 
radiation equality, when perturbations begin to grow mainly at the expense of cold 
dark matter like WIMPs (see, e.g. Ref. [3T]). 
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